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We consider the problem of an harmonic oscillator coupled to a scalar field in the framework of 
recently introduced dressed coordinates. We compute all the probabilities associated with the decay 
process of an excited level of the oscillator. Instead of doing direct quantum mechanical calculations 
we establish some sum rules from which we infer the probabilities associated to the different decay 
processes of the oscillator. Thus, the sum rules allows to show that the transition probabilities 
between excited levels follow a binomial distribution. 

PACS Number(s): 03.65.Ca, 32.80.Pj 

I. INTRODUCTION 

Recently, in analogy with the renormalized fields in quantum field theory, the concept of dressed coordinates has 
been introduced [1-3]. This concept was introduced in the context of an harmonic oscillator interacting linearly with 
a massless scalar field and allowed the construction of dressed states, with the physically correct property of stability 
of the oscillator ground state in the absence of field quanta (the quantum vacuum). This concept also showed the 
technical advantage of allowing an exact nonperturbative treatment of the oscillator radiation processes. 

Indeed, in spite of the system being linear, if we work with the usual coordinates it is not possible an analytical 
computation of the probabilities associated with the oscillator radiation processes. For example, the exact calculation 
of the simplest process, the probability for the oscillator to decay spontaneously from the first excited level to the 
ground state is not possible. However in the framework of dressed coordinates, it has been possible to solve this 
problem exactly and for small coupling constant, it agrees with the long time know result: e~ r * [1]. Furthermore, 
when applied to a confined atom (approximated by the oscillator) in a spherical cavity of sufficiently small radius [2], 
the method accounted for, the experimentally observed, inhibition of the spontaneous decaying processes [4,5]. Also, 
in Refs. [6,7] the case of the nonlinear interaction between the oscillator and the field modes was treated. 

Nevertheless, in all previous works [1-3] it was only computed the probability associated with the simplest process 
just above described. The aim of this paper is to fill this gap. In doing this, we introduce a calculational method 
that we believe could be extended to other physical situations. Instead of performing direct computations of the 
quantum mechanical formulae we establish some identities, that we named as sum rules, and from them we infer the 
probabilities associated with the different oscillator radiation processes. 

The paper is organized as follows. In section II we review the concept of dressed coordinates and dressed states. 
Section III is devoted to the direct computation of the probabilities associated with simplest decaying processes of 
the oscillator. In section IV we introduce the sum rules. Finally, in section V we give our concluding remarks. Along 
this paper we use natural units c = fi = 1. 

II. DRESSED COORDINATES AND DRESSED STATES 

In this section, in order to make this paper self contained, we review the concept of dressed coordinates and dressed 
states as defined in Refs. [1-3]. For this end we consider the system composed by an harmonic oscillator linearly 
coupled to a massless scalar field. The Hamiltonian for this system is given by [1] 
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where qo is the oscillator coordinate, q k are the field modes; k = 1,2,...; Wfe = 2n/L, c k = i]uj k , i) = \ / 2gAuj 1 
Auj = ujk+i — u>k = 2tt/L, g is a frequency dimensional coupling constant and L the diameter of the sphere in which 
we confine the oscillator-field system. In Eq. (1) the limit N — > oo is understood. The last term in Eq. (1) guarantees 
the positiveness of the Hamiltonian and it can be seen as a renormalization of the oscillator frequency [8,9]. 

To understand the necessity of introducing dressed coordinates and dressed states let us consider the following 
problem. Imagine that the oscillator of frequency luq in Eq. (1) represents an atom and the other oscillators represent 
the modes of the electromagnetic field. If there are no interactions among them the free Hamiltonian, c k = 0, has the 
following eigenfunctions, 
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The physical meaning of ip no n 1 ...n N (q) m this case is clear, it represents the atom in its n -th excited level and n k 
photons of frequencies u> k . Now, consider the state tp no o...o{q) : the excited atom in the quantum vacuum. We know 
from experience that any excited level of the atom is unstable. The explanation of this fact is that the atom is not 
isolated from interacting with the quantum electromagnetic field. This interaction in our toy model is given by the 
linear coupling of qo with qk- Obviously, when we take into account this interaction any state of the type t/Vi o...o(<z) 
is rendered unstable. But, there is a problem, the state ipoo...o{q), that represents the atom in its ground state and 
no photons, is also unstable contradicting the experimental fact of the stability of the atom ground state. What is 
wrong? The first thing that cames in our mind is to think that the model given by Eq. (1) is wrong. Certainly, we 
know that the correct theory to describe this physical system is quantum electrodynamics. On the other hand such 
a description could be extremely complicated. If we aim to maintain the model as simple as possible and still insist 
in describing it by the Hamiltonian given in Eq. (1) what we can do in order to take into account the stability of 
the atom ground state? The answer lies in the spirit of the renormalization program in quantum field theory: the 
coordinates q^ that appear in the Hamiltonian arc not the physical ones, they are bare coordinates. We introduce 
dressed (or rcnormalized) coordinates, q' and q' k , respectively for the dressed atom and the dressed photons. We 
define these coordinates as the physically meaningful ones. In terms of these coordinates we define the dressed states 
by 
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where the subscript d means dressed state. The dressed states given by Eq. (3) are defined as the physically measurable 
states and describe in general, the physical atom in the no-th excited level and n k physical photons of frequency w k - 
Obviously, in the limit in which the coupling constant c k vanishes the renormalizcd coordinates q'^ must approach the 
bare coordinates q^. Now, in order to relate the bare and dressed coordinates we have to use the physical requirement 
of stability of the dressed ground state. The dressed ground state will be stable only and only if it is defined as the 
eigenfunction of the interacting Hamiltonian given by Eq. (1). Also the dressed ground state must be the one of 
minimum energy, that is, it must be defined as being identical (or proportional) to the ground state eigenfunction 
of the interacting Hamiltonian. From this definition, one can construct the dressed coordinates in terms of the bare 
ones. Then, the first step in order to obtain the dressed coordinates is to solve the eigenfunctions of the Hamiltonian 
given in Eq. (1). This bilinear Hamiltonian can be diagonalized by introduzing normal coordinates and momenta Q r 
and P r , 
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where {t 7 ^} is an orthonormal matrix whose elements are given by [10], 
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with fl r being the normal frequencies corresponding to the collective modes of the coupled system and given as 
solutions of the equation 
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In terms of normal coordinates and momenta the Hamiltonian given by Eq. (1) reads as 
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then, the eigenfunctions of the Hamiltonian can be written as 

<f>n n 1 ...n N (Q) = (Q\no,Tll, ...,n N ) c 



N 

n 

r=0 . 



1/4 



H„ r (^/Q r Q r )e 2 s 



(8) 




where the subscript c means collective state. Now, using the definition of the dressed coordinates: ipoo...o(q') oc 
0oo. ..a(Q), and using Eqs. (3) and (8) we get e^J^M"'^) = e~2 £™=o &rQl ; f rom which the dressed coordinates 
are obtained 
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We have to point out that the dressed coordinates here introduced are not simply a change of variables, they are 
new coordinates in its own right and are introduced by physical consistence requirement of the model. Also we have 
to stress that our dressed coordinates are not the same as the ones employed in other references, as for example in 
[11] and references therein, where the authors called dressed coordinates the collective normal ones. Also our dressed 
states are different from the ones defined in Refs. [12-16], where the authors called dressed states the states obtained 
by diagonalizing the Hamiltonian of a two level atom coupled to an finite number of electromagnetic field modes. 

Before leaving this section we would like to remark that a similar model to the one given by Eq. (1) has been used 
repeatedly from time to time as a simplified model to describe the quantum Brownian motion [17-20], the decoherence 
problem and other related problems [21]. Since the focus of these works is other the one we have posed here, it is 
not necessary the introduction of the dressed coordinates. On the other hand in quantum optics, the model obtained 
from Eq. (1) by considering the rotating wave approximation has been extensively used [15,22-24]. For example, a 
common situation that is described with this model is the problem of a cavity mode (described by the oscillator of 
frequency wq) coupled to the external modes of the cavity. In this case, the introduction of the dressed coordinates is 
unnecessary. The reason for this is that, in the rotating wave approximation the linear interaction contains only terms 
like a a,k and a k ao and since ajo,fc|0,0, ...,0) = and a^,do|0,0, ...,0) = it is guaranteed automatically that ground 
state of the non interacting Hamiltonian is also eigenstate of the interacting Hamiltonian, assuring in this way the 
stability of the state with no photons both inside and outside of the cavity. Of course if we use the model given by 
Eq. (1), with no rotating wave approximation, to describe the mentioned physical situation it will be necessary the 
introduction of the renormalizcd coordinates, otherwise the state of no photons both inside and outside the cavity will 
evolve to an state of non zero photons in contradiction with experiment. The advantage of our approach in treating 
this problem will be the avoiding of the rotating wave approximation and also, from a technical point of view the 
calculations will be greatly simplified. In particular, this will be seen in next sections where we compute easily the 
probabilities associated with the different oscillator radiation processes. 



III. THE DECAY PROCESSES 



We are interested mainly in the computation of the probabilities associated with the different radiation processes 
of an excited state of the oscillator. Thus, wishing to maintain the reasoning as general as it is possible we show 
the necessary steps to compute the probability amplitude associated with the most general transition: Let the initial 
state of the system, at t = 0, given by |n , n\, n^d, then, we ask what is the probability amplitude of finding it at 
time t in the state \m n ,mi 7 ...mjv),i? Such probability amplitude, which we denote as being A™®™ 1 "^ N (t), is given 

by ' 
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Above, the wave functions <fii i 1 ...i N (Q) and ip no n 1 ...n N (q') are normalized in Q coordinates. The wave function 
<l>loli..-l N (Q) as given in Eq. (8) is already normalized. On the other hand the dressed wave function ip na n 1 ...n N {q l ) as 
given in Eq. (3) is normalized in the q' coordinates but not in Q coordinates. It is easy to show that if the dressed 
ground state ipoo...o(q') is normalized in Q coordinates then, automatically all the dressed states, given by Eq. (3), 
are also normalized. Therefore, assuming that it is the case, we can replace Eqs. (3) and (8) in Eq. (11) and using 
Eq. (9) we get 
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To compute the above integral it will be useful the following identity [25], 
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Next, we compute the probability amplitude of an initial state |00...n jU ...0)d in t = to be found after a time t in 
the state |00...TO„...0)<j. From Eq. (10) we get 
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And after using Eq. (13) in (12) we obtain 



f=0 



ploll...lN 

t OO...n u ...O — 



E 



N\s N 



^+s 1 +...s N =n li 



Substituting Eq. (15) in Eq. (14) we get 
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By setting /j, = v = in Eq. (16) we obtain the probability amplitude for the oscillator to remain at time t in the 
n-th excited state, thus, it reads 

A n Z::° (t) = [Mi)]" , (is) 

where we have discarded the phase factor e ~ ltEaaa -- because it does not contribute for the associated probability. Also, 
by setting /i = and v = k we can obtain the probability amplitude of the particle to decay from the n-th excited 
level to its ground state by emission of n field quanta of frequencies (or energy in ft = 1 units) u>k, 

^°oo.T° = [fok(t)] n . (19) 

In the continuum limit, L — ► oo, the quantity foo(t)j obtained from Eq. (17) by setting /i — v — 0, has been computed 
in Ref. [1]. It reads as 
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For g/wo <C 1, that corresponds to weak coupling, the integrand in above equation is sharply peaked around Q = cj 
and in this case we can obtain easily foo(t) = e~ lLJot ~ rt / 2 , with T = ng. Replacing this result in Eq. (18) and taking 
the square modulus we obtain for the probability that the particle oscillator remains in the n-th excited level, the old 
know result: e~ nTt [26]. 



IV. SUM RULES 



It is clear that if the oscillator is initially in its n-th excited level it can decay of many different ways from that 
described in Eq. (19). For example it can decay to its ground state by emission of n field quanta of different frequencies 
or it can also decay to other lower excited states by emission of a number of field quanta less than n. The probability 
amplitudes related to these processes can all be computed by using Eqs. (10) and (12). The task to be made can be 
very hard because it would be necessary to compute the integral given in Eq. (12) that involves products of more 
than two Hermite polynomials, as it can be noted by substituting Eq. (13) in (12). Therefore, we will avoid this 
complication following an alternative way. If the particle oscillator is at time t = in the n-th excited state we expect 
that at time t: it remains excited, emit a field quantum of frequency u>ki and go to the (n — l)-th excited level, emit 
two field quanta of frequencies u)^ , ujk 2 an d go to the (n — 2)-th excited level, emit three field quanta of frequencies 
LUk ± , Uk 2 , t0k 3 and go to the (n — 3)-th excited level and so on. We denote the probability amplitudes related with 

these processes respectively by ^oo. 1 .'o* 1 (*)> •^ioo.^o'' 11 * 2 ^o.^o'" 1 ^ ^ (*)> or in general A^^l" 11 " 2 "' 1 "' (t), 
i = 1, 2, ...n. The corresponding probabilities are denoted by "P^o 1 ^ kllk2 " Aki (f) ) i = i ) 2...n. Now we will compute 
all these quantities from the knowledge of the probability of the particle to remain excited in the n-th excited level, 
"P™oo V.o (*) whose probability amplitude is given by Eq. (18). For this end we use the identity 

Ko°o°7.o° W + E V no^o kl (*) + E ^oo. 2) o lfel lk2 (*) + - + E Klo 1 . 1 .o 2 - lkn (*) = 1 , (21) 
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that expresses, the sum of the probabilities of all possibilities is equal to one. Starting from Eq. (18) we have to be 
able to get an identity similar to Eq. (21) and from such expression identify the respective probabilities associated 
with all the other possible processes. 

From Eq. (18) and using Eq. (17) with fi = v = we can write ^"oo.'.'.oC*) as 
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By adding and subtracting 1 on the right hand side of the Eq. (22) and using the relation given by Eq. (23), we 
obtain 
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From the orthonormality of the matrix elements {t 7 ^} we have X^^a* = t^tf) + ^ k t k t k — 8 rs and using it in Eq. 
(24) we get 
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where in the second line we have used the symmetry of the expression under index permutations. In terms of /oo(i) 
and /ofe, that can be obtained from Eq. (17), we can write Eq. (25) as 
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an identity of the type we are looking for, compare with Eq. (21). This identity is what we call in the present paper 
as sum rules. Other similar identities can be established for other related processes. Our sum rules are very different 
from the quantum mechanical sum rules [28] where some identities are established by the only use of the algebra 
between canonically conjugated variables. In the present case we used the same name since our sum rules can be 
established by the only use of the algebra of the matrix elements that diagonalize the Hamiltonian. 

From Eq. (26) we can identify the respective probabilities associated with all the other possible processes, thus, by 
comparing it with the Eq. (21) we obtain 
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from which we can also obtain the corresponding probability amplitudes which are given as 
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As a check, that expression above is a valid expression, we set i — n and k\ = k 2 = ... = k n = k in Eq. (28) and we 
obtain Eq. (19) which is the probability amplitude to the emission of n field quanta of frequencies uj k as it must be. 

From Eq. (27) we can also obtain the probability of the dressed oscillator to decay, at the time t, from the n-th to 
the m-th excited level by emission of i = (n — m) field quanta of arbitrary frequencies. We denote this quantity by 
V n ^m{t)- For this end, in Eq. (27), we sum up over all possible values of hi, fc 2 ,...,fci and using the identity 



J2\Mt)\ 2 = l-\foo(tf 



(29) 



that is obtained from the orthogonality property of the matrix elements we get 

Vn^m(t) = Ul r. (\foo(t)\T (l - |/oo(t)| 2 ) , m = 0, 1, ...,n; (30) 

m!(n — mj! 

where we can note that V n ^m(t) is given by a binomial distribution [29]. 

To our knowledge no similar result to Eq. (27) or (28) has been obtained previously. On the other hand, for weak 
coupling |/oo(i)| 2 = e ~ Tt i we can write Eq. (30) as 

7>„W - ,, "' „ e- rort (l - e- Tt T- m , m = 0, 1, n; (31) 
m!(n — m)! 

a result similar to the one obtained in the early days of quantum mechanics [27] (in the notation of this reference 
i = n — in) by a method similar to the Einstein derivation of the black-body radiation formula. 

V. CONCLUDING REMARKS 

The probability amplitudes for other processes, given by Eq. (10), can be obtained by using the crossing relation 
present in that equation. We can see from Eq. (10) that the probability amplitude of an initial state \no,ni, ...,njv)d 
to be found at time t in the state \mo, mi, mjv)d is the same as the probability amplitude of the initial state 
|m ,TOi, ...,miv)d to be found at time t in the state \n ,ni, ...,njv)d- In particular this means that the probability 
amplitude of emission of a field quanta is the same as the probability amplitude of absorption of the field quanta. We 
have just computed, Eq. (28), the probability amplitudes related to the emission of field quanta. Using the crossing 
symmetry mentioned we can compute all the probability amplitudes related to the absorption of field quanta. The 
crossing relation mentioned is different from the field theoretical crossing relations, where the probabilities are equal 
and not the probability amplitudes as it is in the present case. The reason for the occurrence of this in our present 
model is because the wave function of an harmonic oscillator is real. 

Thus, we have shown that any probability amplitude associated to the radiation processes of an harmonic oscillator 
in interaction with a massless scalar field can be given in terms only of two quantities: the probability amplitude of the 
particle oscillator to remain in the first excited level foo(t) and the probability amplitude of the particle oscillator to 
decay spontaneously by emission of a field quantum of frequency u>k, fok(t)- We have also shown that the probability 
for the spontaneous decay of the oscillator from one excited level to a lowest one, by emission of arbitrary field quanta, 
is given by a binomial distribution. 

Finally we have to stress that in computing the probabilities associated to the above mentioned processes, the 
calculations are greatly simplified. This can be noted particularly in computing the integral given by Eq. (12). If no 
dressed coordinates were introduced, the integral that would appear instead will not contain the exponential factor 
e~^r=o il rW r ^ but other more complicated term that will prevent us from using directly the orthogonality properties 
of the Hermite polynomials. In this way no exact calculations will be possible. Then, we believe that the use of 
the dressed coordinates concept will greatly simplify the study of early works, where extensive use has been made of 
the model with Hamiltonian given by Eq. (1) to model different physical situations, such as the quantum Brownian 
motion, decoherence and other related problems in quantum optics. The study of these problems, in the framework 
of dressed coordinates, are under study and will be reported elsewhere. 
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